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13. ABSTRACT 

The report presents results of continued resear.ch aimed at improving the predictability of 
.off-nominal internal ballistics performance of solid propellant rocket motors (SRMs) 
including thrust imbalance between tvo SRMs firing in parallel. The previous design 
analysis and Monte Carlo computer programs developed by the authors for use in such perform- 
ance analyses are extensively revised to permit easier use and substantial reduction in 
computer operating time. An earlier analysis of internal ballistic effects of circular- 
perforated propellant grain deformation is verified by additional comparisons with test 
results. The potential effects of nozzle throat erosion on internal ballistic performance 
I are studied and a new propellant burning rate law postulated. The proonllant burning rate 
j model when coupled with the grain deformation model permits an exceller.t match between 
theoretical results and test data for the Titan II IC, TU455.02, and the first Space Shuttle 
SRM (OM-1). Analysis of star grain deformation using an experimental model and a finite 
element model shows the star grain deformation effects for the Space Shuttle to be small in 
comparison to those of the circular-perforated grain. An alternative technique is developud 
for predicting thrust imbalance without recourse to the Monte Carlo computer program. A 
scaling relationship used to relate theoretical results to test results may be applied to 
the alternative technique of predicting thrust imbalance or to the Monte Carlo evaluation. 
Extended investigation into the effect of strain rate on propellant burning rate leads to 
the conclusion that the thermoelastic effect is generally negligible for both steadily 
increasing pressure loads and oscillatory loads. The capability to calculate the inflnenro 
of oscil laying pressure loads on transient burning rates making use of the Zeldovich and 
Novozhilov transient burning rate model demonstrated during the course of this investi.iaiioi: 
,may prove useful in combustion instability investigations. 
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NOMENCLATURE 


English 

Symbol Definition 


a Propellant burning rate coefficient, 

a^ Time-temperature shift factor. 

A Cross-sectional area. 

b Average web- thickness of main propellant grain, 

c Specific heat at constant volume. 


c. Constant giving ratio of peak amplitude of cyclic 

pressure loading to average pressure. Also, 
constant in the Summerfleld burning rate law. 

C 2 Constant in the Summerfleld burning rate law. 

Cj,c^ Constants in modified burning rate law. 

Cj^.C^ Constants in time-temperature shift factor (Eq. V- 

C„ Ideal thrust coefficient. 

F 

o 

CV Coefficient of variation. 


D The.rmoelastlc energy dissipation. 

E,E* Modulus of elasticity and real part of the complex 
tensile modulus, respectively. 




Modulus of elasticity under long-time steady 
loading conditions and at the glass transition 
temperature, respectively. 


F Thrums t. 


G Mass flow rate per unit area. 

G' Real part of the complex shear modulus. 

K 2 Two-sided statistical tolerance factor. 

L,L Flow length and reference flow length, respec- 

tively, in modified flame height burning rate law. 


ix 


Units Used 
in/sec-psi” 

in2 

in. 

in-lbf/lbm“R 
,lbf-sec/in^ 

lbfl/3_sec/inS/3 
sec, (sec^/slug) 

. /K 

®R/sec 

lbf/in2 

lbf/in2 

Ibf 

slugs/sec-in^ 

Ibf/in^ 

'in. 
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NOMENCLATURE (Continued) 


English 


Definition 


Units Ul'#- 


Mach number. 


Burning rate exponent. 


Pressure. 


Critical pressure in modified flame. 
Burning rate. 

burning rate at zero velocity. 


In/sec 

In/sec 


Reynolds number. 


Standard deviation of a sample. 


Time. 


Absolute temperature 


’R or *K 


Glass transition temperature. 

Velocity parallel to the burning surface. 


ft/sec 


Difference in web thickness at ends of L- 
(Ref. 5). ® 

Value of general statistically distributed 
variable. 


Distance burned. 


Greek 

Symbol 


Instantaneous difference between web thickness 
at head and nozzle ends of controlling grain 
length, Including any length associated with 
X.j.^ and 0Q (Ref. 5). 


Linear coefficient of thermal expansion. 


Shear loss and tensile lo.^s modulus, respectively. 


X 



NOMENCLATURE (Continued) 


6r«ek 

Symbol Definition Units Used 

$ Constant exponent used In rolsxatlon modulus 
expression. 

Y Ratio of specific heats. 

A Change or difference in quantity. units vary 

t Strain. 

Cp Thrust loss coefficient. 

X Thermal conductivity. 

V Poisson's ratio, 

p Density. 

0 The standard deviation of a statistically 

distributed variable; i.e., the square root 

of the second moment about its mean value. 

Also peak value of oscillatory normal stress. 

t Web thickness. Also peak value of the 

oscillatory normal stress. 

Web thickness. 

T Time constant used in relaxation modulus 

^ expression. 

(I) Cyclic frequency. 

Subscripts Definition 

a Ambient. Also burning rate coefficient.. 

av Average or mean value. 

e Exit plane of notsle. 

h Head end of grain. 

n Nosale end of grain at position of maximum flow Hach 

number. 


in**lbf/ in-sec *R 

Ibm/in^ 

units vary 
Ibf/ln^ 

In-lbf/ln^ 

la. 

sec. 

H* 
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NOMENCLATURE (Concluded) 


Subscripts Definition 

N Value for new motor, 

o Stagnation. 

S Space Shuttle. 

SA Adjusted value for Space Shuttle, 

ss Steady state, 

t Transient . 

T Titan IIIC. 

Superscripts 

• Time rate of change of a quantity. 

' Experimental data. 

* Value of a variable where M»l. 


Average or mean value. 



I. INTRODUCTION AND SUMMARY 


This report presents the results of research performed at Auburn 
University during the period November 19, 1976, to November 18, 1977, 
under Modification No. 19 to the Cooperative Agreemc.it, dated February 
11, 1969, between NASA Marshall Space Flight Center (MSFC) and Auburn 
University. The principal r ijectlve of the research is to further 
develop techniques for theoretical assessment of solid rocket motor 
(SRM) internal ballistic performance to include statistical investigation 
of thrust Imbalance of pairs of SRMs firing in parallel as on Che booster 
stage of the Space Shuttle. 

The theoretical thrust imbalance of SRM pairs has been previously 
Investigated statistically using a Monte Carlo technique (Ref. 1-4). 

The results of this investigation include a computer program which selects 
sets of significant variables on a probability basis and calculates the 
performance characteristics for a large number of motor pairs using a 
mathematical modei^ of the internal ballistics. This model is based on an 
earlier basic computer program (Refs. 5 and 6) which permits rapid assess- 
ment of internal ballistic performance using simplifying assumptions which 
have been demonstrated to give realistic results. The basic program is 
known as the design analysis program. A complete listing of the design 
analysis program is given in Ref. 3. A few minor modifications are 
identified in Ref. 2. Reference 2 also lists the complete program with 
a special option which permits assessment of the effects of circular- 
perforated grain deformation on internal b.iiilstics. As adopted into the 
Monte Carlo program, the basic design analysis program does not Include 
the grain deformation analysis but docs Include an analysis which accounts 
for the effects of propellant grain out-of-roundness , eccentricity and/or 
nonaxlsyinmetric grain temperature distribution. A listing of the complete 
Monte Carlo program is also given in Ref. 3 with a few minor improvements 
identified in Ref. 2. 

The research which culminated in the development of the design analysis 
and Monte Carlo programs revealed a number of areas in which the programs 
may be improved. Also, an alternative scheme for predicting thrust im- 
balance has been suggested, and the methods of assessing grain deformation 
and its ballistic effects need further verification. The extended research 
discussed in this report deal.s in part witli such additional aspects. 

During development of the Monte Carlo program the effects of strain 
rate on the performance of SRMs were examined (Refs. 2-3). The investiga- 
tion disclosed that modification of the solid phase grain temperature 
distribution and hence propellant burning rate due to thermoelastic coupling 
may occur when strain rates are high; i.c., during highly transient con- 
ditions of operation. Although an assessment of single SRMs under ignition- 
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type pressure loads showed the burning rate effect la smell* additional 
research was Indicated with respect to effects of physical properties* 

Also* the thermoelastic effect needed to be evaluated under cyclic loading 
conditions such as those produced by oscillatory burning and mechanically 
Induced vibrations. The extended research also embraces Investigation of 
these areas. 

The alternative scheme for predicting thrust isibalance that resulted 
from the extended research makes use of a demonstrated linear relation- 
ship between thrust Imbalance at various times and the corresponding slopes 
of the nominal thrust-time trace for both theoretical (Monte Carlo) results 
and test data. A method Is presented whereby the thrust imbalance of a new 
SRM design may be predicted from the linear relationship without recourse 
to a Monte Carlo analysis. The accuracy of the prediction Is expected to 
be best during the talloff portion of the trace. The method la restricted 
to SBMa employing similar manufacturing and processing controls to those of 
the Titan llIC SRMa which constitute the baseline analysis. An exception 
Is propellant burning rate variability which may be different for the new 
design since It Is accounted for In the relationship. A scaling relation- 
ship Is used to relate theoretical results to test results. This scaling 
relationship may be applied either to the alternative technique of pre- 
dicting thrust imbalance or as an adjustment to a Monte Carlo analysla of 
thrust Imbalancn. 

Improvements to the basic design and Monte Carlo computer programs are 
made In the following ateas: 1) Modification of the input procedure In the 

programs to allow the use of NA^(F.I.IST for nonstatistlcal Input variables. 

This simplifies the mechanics of card input for both programs and reduces 
the amount of Input required for multiple configurations In the design 
analysla program. 2) Incorporation of an output option In the design 
analysis program which allows selection of a limited number of output var- 
iables greatly reducing the number of printed sheets and computer execution 
..Ime. 3) Reduction in the volume of printed output In the Monte Carlo 
program by changes in the present format structure. This change causes all 
Input data for each SRM to be written on a single page. A) Modification 
of program logic and statement structure in both programs to decrease 
computational time and hence further increase program efficiency. The net 
result of all the changes Is a redaction of typical computational times In 
the basic design program by approximately AO percent for the limited out- 
put option and 20 percent for the full output , Monte Carlo program opera- 
tional times are reduced by approximately 20 percent. 

Attempts to further verify the influence of grain deformation on 
internal ballistics were Impeded by the anomalous performance of the motor 
selected for the Investigation. Efforts to explain the Irregular behavior 
Include analysis of nosile throat behavior and propellant eroaive burning 
characterlatlca. A new model of eroaive burning, termed the modified flame 
height model* Is developed which when combined with the grain deformation 
model permits excellent representation of the performance of the Caator 
TX35A-5 and TU-A55.02 SRMs. Tho model uaoa two constants which were selscted 
for each motor to match test performance. The same constants used for the 
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TU-4S5.02 also gave a good match when used to predict the performance of 
the first Space Shuttle SRM (DM-1) which has a very similar, although not 
identical, propellant iormulation. The new burning rate model was in- 
corporated as an option into the same version of the design analysis pro- 
gram that contains the grain deformation analysis along with the other 
design analysis program Improvements outlined above. Although this research 
provides additional evidence of the validity of the grain deformation 
analjFsls for circular-perforated grains, we feel that a sufficiently firm 
basis has not yet b*den established to warrant incorporating the analysis 
into the Monte Carlo program. 

With the objective of finding a simplified approach to Including the 
effects of star-type grain deformation into the program, an experimental 
subscale inert model of the Space Shuttle star grain segment was constructed. 
Perimeter changes measured under hydraulic pressure were approximately 
0.2 percent at 400 psla. This was confirmed by a finite element analysis. 

The potential amount of burning perimeter change involved was Judged to be 
small enough in comparison to circular-perforated changes (4.0 percent at 
800 psi) to be considered negligible in the Space Shuttle. However, for 
grains with thicker webs and shorter star points the deformation is expected 
to be more substantial. 

Additional investigation into the effect of propellant strain rate on 
propellant burning rate leads to the conclusion that for all practical 
purposes the effect is negligible for both steadily increasing pressure loads 
and oscillatory loads. The investigation includes incorporation of a relax- 
ation modulus into the thermoplastic analysis and examination of the effect 
of Poisson’s ratio, thermal expansion and thermal conductivity. A brief 
analysis of the ef*’ects of mechanical vibration on burning rate and internal 
ballistic performance as caused by the thormomechanlcal energy dissipation 
associated with viscoelastic materials is also Included. An important by- 
product of this research is the demonstrated capability to cr.lculace the 
influence of oscillating pressure loads (with or without thermoelastic 
coupling) on transient burning rates making use of the Zeldovlch .and 
Novo*hllov (Z-N) transient burning rate model (Refs. 2 and 7). Thi.s analysis 
may prove useful for combustion instability studies Irrespective of the 
thermoelastic effect. Likewise, an Internal ballistic model of the Ignition 
transients in SRMs also including the .Z-N burning rate model developed in a 
separate effort (Ref, 8) under the general guidance of this project may prove 
of some utility. 

As well as a general discussion of the changes that have been made to 
the basic computer programs, the text contains a detailed description of 
the new input formats for the design analysis and Monte Carlo prograti's 
and for the design analysis program with grain deformation and burning rate 
model options. Sample outputs are also presented lor the three programs. 



II. THRUST IMBAUNCE PREDICTIONS 


In Ref. 2, a statistical relationship is shown to exist between the 
thrust imbalance at various times and the corresponding elope of the 
nominal thrust-time trace for 30 theoretical Titan IIIC SRH pairs obtained 
from a Monte Carlo evaluation. The purpose of the present analysis Is to 
refine and extend the investigations and conclusions of Ref. 2. More 
specifically, the new analysis ostabIi.shes procedures for predicting the 
thrust imbalance trace for a new SRM pair directly from the slope of the 
nominal thrust-time trace without a Monte Carlo evaluation. Also, a 
scaling law is Uovclopcd based on differences between the Monte Carlo 
evaluation and known lost results for the Titan 1110. 

Linear Approx! im i t ton to Monte Carlo Imb al.'tnce Equations 

A question was raised in Kef. 2 concerning the randomness of selection 
of the time .slices on which the thrust imbalance relationship was based. 

A rigorous statistical analysis of the lelationshlp between the thrust 
imbalance and thu thrust-time slope requires that a random selection of 
slopes be made, l.e., with no prior knowledge of thrust Imbslsncs values. 
This necessitates a special search procedure since for a given value of 
thrust-time slope, several corresponding thrust imbalance values (at various 
times) may be found. The task becomes exceedingly complicated with rela- 
tively low values of thrust-time slope, especially when examining erratic 
tost data. , 

An alternative approach is to divide the thrust trace into a large 
number of equl-spaccd time slices. While tills does not represent s truly 
random selection of thrust-time slopes, it can bo argued that the slopes 
are chosen with no prior knowledge of the thrust Imhalancs values. Although 
not completely rigorous, this procedure provides an accurate end consistent 
method of analysing and predicting thrust Imbalance, as will be demons tr. tied. 

The alternative approach is applied to a set of 30 theoretical Titan 
me SRM pairs. Time slices of motor pair Imbalance traces and individual 
motor thrust traces are taken at 0.3 soc. Increments throughout burn time, 
including tailoff. Average values of absolute thrust imbalance for the 30 
pairs at each time slice are plotted versus the absolute values of ths 
average ihrust-tlme slopes for the bO motor.'*. This is shown in Fig. Il-ls 
alorg with a linear least-squares approximation, 5172. 

A correlation coefficient, which indicates the accuracy of ths linear 
approximation, is determined to bo 0.983. These values agree very well 
with those presented in Ref. 2, 

it .should he noted here that relatively large values of thruet-tlme 
slope and thnist imhal.nu'o .u'c oxpcrl«'nccd only tlurlng the tailoff period. 
Thus the slope of the linuar approxlnat ion Is determined principally by 
the tv^lloff characteristics, ami tho m.wlmum value of thrust imbslsncs 
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ABSOLUTE VALUE OF MEAN TIME RATE OF CHANGE OF THRUST (LBF)xlO” 


a. Thrust Imbalance vs. thrust-time slope 



b. Standard deviation of thrust imbalance vs. thrust-time slope 


Figure II-l. Correlations determined from a set of 30 theort't- 
Ical Titan IIIC SRM pairs obtained from a Monte 
Carlo evaluation. 
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occurs during the tailoff period. Thl.*? observation tends support the 
validity of using equl-spaced time slices in place of random values of 
thrust-time slope since the Inherent, additional quantity of lower valued 
thrust-time slopes, obtained in the equl-spaced time slice method, only 
confirms the thrust lnd>alance Intercept, and has little effect on the 
slope of the established relationship between thrust Imbalance and thrust- 
.me slope. 

It is desired to compare this analysis of 30 theoretical Titan IIIC 
SRM pairs to an analysis of 30 theoretical Space Shuttle SRM pairs. After 
obtaining the latter SRM pairs from a Monte Carlo evaluation, the pro*- 
cedure described above is used to determine the thrust imbalance versus 
thrust-time slope correlations presented in Figure Il-2a, with a linear 
least-squares aonroxlmatlon, “ 0.232 1^1 + 1504, and a correlation 

coefficient of 0.943. It is immediately obvious that the linear approxi- 
mations for the two analyses are In considerable disagreement, which might 
well be expected when comparing two entirely different families of motors. 

It should be noted that the Space Shuttle SRM configuration used in the 
analysis differs somewhat from that used in Ref. 3 and includes changes 
that reduce tailoff thrust Imbalance. The earlier configuration would 
exhibit higher thrust Imbalances at correspondingly higher nominal thrust- 
time slopes. These sets of values would presumably lie close to an 
extension of the regression line on Fig. II-2a. 

In examining the cause of the apparent discrepancy between the Titan 
and Space Shuttle SRM evaluations, It is advantageous to consider the 
sources of the data used; namely, the Monte Carlo evaluations, which use 
statistical data for determining variations of certain parameters. An 
examination of the statistical variations (standard deviations) of these 
parameters and their relative influences reveals that the variation of the 
burning rate coefficient may represent the most significant difference 
between the two analyses. 

The other statistically varvlng parameters were found to have small 
influences on the thrust Imbalance relationships because of similar 
coefficients of variation (standard deviation divided by mean value) for 
these two SRMs. Other SRMs might exhibit different variations of parameters 
and the effects could become influential . 

The coefficient of variation of the Titan burning rate coefficient 
used in the analysis is 


CV_ - /a_ - 0.0003428/0.06653 - 0.005153. 

T *T ^ 


(II-l) 





■ , -^“1 


♦ ♦ ♦♦ ♦ ♦ 


AF =0.232 
' 'av 


+ 1504 


Correlation Coefficient » 0.943 


i20 ikO 200 2^10 260 3*20 

ABSOLUTE VALUE OF MEAN TIME RATE OF CHANGE OF THRUST (LBF)xlO~3 

a. Thrust Imbalance vs. thrust-time slope. 


52- S|^j.|-0.144 iF^^I + 1614 

Correlation Coefficient » 0.935 

o 

0 '*<* «b t^o ilo 3bo alio 55o aio 

ABSOLUTE VALUE OF MEAN TIME RATE OF CHANGE OF THRUST (LBF)xl0“^ 

». Standard deviation of thrust Imbalance vs. thrust-time slope. 


Figure II-2. Correlations determined from a set of 30 theoretical 
Space Shuttle SRM pairs obtained from a Monte Carlo 
evaluation. 
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The coefficient of variation of the Space Shuttle burning rate coefficient 
used in the analysis is 





y, 






CV_ - o /a. - 0.000008335/0.006915 - 0.001205. (11-2) 

b ag b 

And, since the magnitudes of thrust imbalance result primarily from the 
variations of web-action times, not only should the variations in burning 
rate coefficients be examined, but also the web thicknesses. The coef- 
ficient of variation for the Space Shuttle burning rate coefficient may 
now be adjusted to that used for the Titan, including web thickness adjust- 
ment : 


CVgA - CV^/Tg - (36.16 in)(0.005153)/(4C 60 in) - 0.004589. (II-3) 

The adjusted Space Shuttle standard deviation is 


0 • CV-. - (0.006915) (0.004589) - 0.00003173. (II-4) 

®SA ® 

Using this adjusted value of statistical variation for the Space 
Shuttle burning rate coefficient, a Monte Carlo evaluation of 30 theo- 
retical SRM pairs is obtained. Repeating the previously described time 
slice analysis results in the relationship of Fig. ll-3a, with a linear 
approximation, ■ 0*7281 f^^| + 1633, and a correlation coefficient 

of 0.928. The effect of this burning rate adjustment is clearly seen 
when compared, as in Fig. II-4, to the original Space Shuttle and Titan 
analyses. It is noted that the Space Shuttle configuration with adjusted 
burning rate variations closely defines the thrust imbalance versus thrust- 
time slope of the theoretical Titan motors. Thus it is concluded that 
the burning rate coefficient and web thickness are the most significant 
factors influencing the magnitudes of the thrust Imbalance, particularly 
during the talloff period. 

It now appears feasible that a thrust imbalance versus thrust-time 
slope relationship, based on a Monte Carlo evaluation, can be predicted 
for a pair of new motors firing in parallel or a family of new motors, 
without actually performing the Monte Carlo evaluation. However, this 
does require that the new motors have similar manufacturing tolerances 
and statistical variations (similar technology) on the significant motor 
parameters. Using both the Titan and Space Shuttle analyses as a base, 
a thrust Imbalance relationship is determined for the new motors, 

IafI^^ - (0.0714 + 3.282 + 26,694 CV^j t^^) , (II-5) 
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b. Standard deviation of thrust Imbalance vs. thrust-time slope. : 

Figure II-3. Correlations determined from a set of 30 theoretical | 

Space Shuttle SRbt pairs obtained from a Monte Carlo j 

evaluation using Titan burning rate variations. 
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wh«r« CV{{ la Che coefficient of variation of the burning rate coefficient 
tor the nev notors and (In Inches) is the mean web thickness of the 
new tnotora. This relationship is obtained by assuming linear relation- 
ships between the product of CV and t and the slope and Intercept of the 
thrust Imbalanca versus thrust-time slope. 

The discussion to this point has dealt only with thrust Imbalance 
correlations; however, an .Identical procedure has been executed simultan- 
eously to establish correlations for Che standard deviation of the thrust 
Indialauce. Fig. Il-lb represents the standarJ deviation correlation for 
Che Titan theoretical ;>airs; Fig. ll-2b for the Space Shuttle theoretical 
pairs; and Fig. II-3b for the Space Shuttle theoretical pairs using Titan 
variations and web adjustment on burning rate coefficient. These 
analyses are compared in Fig. II-5 to illustrate the effect of the burn- 
ing rate adjustment on the standard deviation. Since the comparison is 
similar to that demonstrated In Fig. II-4 for the thrust loibalance, a 
relationship for the standard deviation of the thrust Imbalance for a 
pair of new motors Is determined, 


*|AF| - (-0.00126 + 2.969 CVj^ Vl^av* % V* 

again without the Monte Carlo evaluation. 

Linear least squares approximations to the actual relationships 
between thrust Imbalance, standard deviation and thrust-time slope have 
been established, but nothing has been said about the applicability of 
these corralatlons. We wish to use these approximations to establish an 
expected thrust Imbalance envelope for the Space Shuttle theoretical SRM 
pairs, as was accomplished In Ref. 2, using the relatlouahlp 

4F . IKj [|4F|*, + (II-7) 


where K 2 Is the two-sided tolerance factor (3.68) determined by Che 
sample else (30 pairs) and the desired probability (0.90) Chat a speci- 
fied percentage (99. 7X) of the thrust imbalanca distribution will be 
Included within the prescribed limits (thrust Imbalance envelope). 

A comparison Is presented In Fig. I 1-6 of the thrust Imbalance 
envelopes using the Monte Carlo evaluations directly and using the linear 
approximations. For conciseness only the positive halves of the envelopes 
are presented; the negative halves are simply mirror Images. Fig. II-6a 
demonstrates a close agreement of the two methods during the talloff 
period. This is a result of the relatively close correlation of the 
linear approximation to the Monte Carlo data during talloff, as seen in 
Pig. 11-2. '' over. Fig. II-6b does not indicate a close agreement of the 
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a. Thrust Imbalance envelopes during tailoff. 


O 



Figure 11-6. Tolerance limits for thrust imbalance of 
30 theoretical Space Shuttle SRM pairs 
(Probability 0.90 for 99. 7S of the dis- 
tribut ion) . 
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two methods during web action time, especially at the end of the head end 
star grain web nctlon time. This discrepancy can be explained since the 
slopes of the liuear approximations were determined principally by the 
talloff characteristics, and the associated greater web action times tend to 
magnify any variations in burning rate at shorter times. The*‘<?fore, the 
slopes of the linear approximations are accurate only during cua talloff 
period. The slope can be adjusted by reducing it in proportion to the 
percentage of web burned at any particular time. Continuous adjustment 
of the slope throughout web action time results in the third thrust im- 
balance envelope represented in Fig. XI-6b. This adjustment does result 
in much closer agreement at the end of star grain web action time, but 
adversely affects most of the remaining envelope. Fince the linear 
approximations are reasonably accurate only for higher values of thrust- 
time slope, it can be argued that the web-time adjustment should only be 
applied to higher magnitudes, such as at the end of the star grain web 
action time. This would, however, present difficulties in transition 
between adjusted and non-adjusted portions of t;he .ivelope. 

Since the linear approximations of the thrust imbalance and the 
standard deviation of the thrust Imbalance are reasonably accurate only 
during talloff, a better approximation of t..e talloff imbalance envelope 
could possibly be achieved by using only the talloff trace to determine 
the linear approximations. 


Because of the somewhat random relationship between thrust imbalance 
and thrust-time slopes during web action time, this portion of the envelope 
should be used only as a first, and generally not conservative, estimate. 

For the talloff period, an alternative method has been devised whereby the 
thrust Imbalance characteristics obtainable from the Monte Carlo program 
can be estimated quite well from the nominal thrust-time trace using the 
linear approximation given by Eqs. II-5 and 6. The equations ace only 
valid, however, for motors having manufacturing and processing variations 
similar to the Titan IHC and Space Shuttle SRM's, except for burning 
rate control, which is accounted for in the equations. 

Comparison with Test Results 

Thrust Inlalance relationships have now been established between 
Monte Carlo evaluations of two different SRM configurations. The question 
now arises as to how these theoretical evaluations compare with actual 
test data. At present, the Space Shuttle SRMs have not flown and no test 
data for motor pairs is available, so we can institute a comparison only for 
the Titan IIIC configuration, for which good paired thrust data is readily 
available for 20 SRM pairs. 

Thrust Imbalance and thrust-time slope are correlated for the 20 
Titan IIIC SRM pairs of test data in an identical manner as performed 
previously for the theoreHcai motor data. Figure Il-7a depicts the 
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b. Standard deviation of thrust Imbalance vs. thrust-tlmc slope. 


Figure II-7. Correlations determined from flight data obtained 
from a set of 20 Titan IIIC SRM pairs. 
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results of this correlation along with a linear least-squares approximation 
IafI » 0.492 If I + 3963, with a correlation coefficient of 0.955. Com- 

parlson of the flight data and theoretical motor correlations reveals sig- 
nificant differences In thrust Imbalance characteristics. Because the 
equations used In the Monte Carlo evaluation do not represent all the complex 
phenomena associated with the beginning of talloff, the resulting theoretical 
traces Initially regress more abruptly than the actual performance traces. 
Therefore, the maximum thrust Imbalance Is higher and occurs earlier In the 
theoretical motor pairs than In the actual pairs which have more rounded 
traces at the beginning of talloff. It Is also notable that whereas the 
Monte Carlo program overestimates the maximum Imbalance based on the present 
time slice analysis. It underestimates the maximum thrust imbalance calcu- 
lated without regard to the time at which that maximum occurs (See Ref. 2). 

Because of the nature of the differences In talloff characteristics 
between the theoretical and actual Titan IIIC motors. It Is reasonable to 
assume that these differences will be exhibited In other configurations 
such as the Space Shuttle SRMs. With this assumption it is now possible, 
using the theoretical Space Shuttle SBM thrust imbalance correlation 
(Fig. Il-2a) , to predict the thrust imbalance correlation of the flight 
article. Using the proportionate changes in slope and intercept between 
the Titan correlations (Figs. II-l and 7) and applying these to the Space 
Shuttle theoretical correlation (Fig. 11-2) , results In the Space Shuttle 
flight prediction as depicted in Fig. II-8 which is = 0.167 

+ 1152. Additionally, it may be noted that a proportionate Increase in 
mflx<mum thrust- time slope Is included in the Space Shuttle prediction, 
since this Increase is also exhibited by the Titan motors during talloff. 

A similar scaling is applied to the standard deviation of the thrust 
Imbalance as Indicated in Figure II-9. The Space ^ Shuttle flight prediction 
of the standard deviation is then “ 0.097 I.F^^j + 1465. 

Using these flight predictions (scaled linear approximations), an 
expected talloff thrust Indjalance envelope Is determined and presented 
in Fig. II- 10a. The validity of the scaled thrust Imbalance envelope can 
be verified by examination of the Titan IIIC theoretical and test data. 

Figure Il-lOb depicts the thrust imbalance envelope computed from: 

1) the Monte Carlo Evaluation, directly; 2) linear approximations of the 
Monte Carlo evaluation; 3) scaled linear approximations (flight predictions); 
and 4) actual test data. It Is noted that the scaled envelope is slightly 
lower than that of the test data. This is because the linear approxima- 
tions of the test data (Fig. II-7) do not accurately represent the maxifjm 
values of thrust Imbalance and standard deviation since these values lie 
off the regression line. Also, In Fig. Il-lOa, the time of peak thrust 
Imbalance differs between the Monte Carlo and the test data analyses, as 
described previously in this report. However, the magnitude of the peak 
imbalance Is fairly accurately represented by the flight prediction 
envelope . 
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Monte Carlo Evaluation (30 pairs) 



a. Space Shuttle thrust imbalance envelopes. 


Monte Carlo Evaluation (30 pairs) 
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The scaling relationship between the theoretical evaluation and 
flight test results may be applied either to the new procedure with no 
Monte Carlo analysis or to a Monte Carlo analysis of a new SRM pair. 

The latter alternative Is recommended when manufacturing and processing 
controls differ substantially from those employed in the Titan IIIC. 
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III. IMPROVEMENT OF PREVIOUS COMPUTER PROGRAMS 

Because of the anticipated use of the design analysis and Monte 
Carlo computer programs In the preliminary design of SRMs, It Is desirable 
at this time to Increase the efficiency of the programs and reduce some 
of the complexities associated with handling of data. This Is accomplished 
through an extensive series of modifications aimed at reducing computer 
execution time, reducing program size, eliminating unnecessary variables, 
statements and computations, reducing the amount and nuinber of cards of 
Input data, allowing more flexibility and easier modification of Input 
data, and reducing the volume of output data. 

Design Analysis Computer Program 

In an attempt to accomplish the aforementioned goals, many modifica- 
tions have been made. The following represents a list of the most signif- 
icant modifications since Ref. 3: 

1. The data card to 'SET INITIAL VALUES OF SELECTED VARIABLES TO ZERO' 

Is removed and this operation replaced by utilizing available 
temporary disk space. 

2. All Input data except NRUNS and tabular data for propellant tem- 
perature distribution are converted to NAMELIST Inputs. This 
facilitates easier data Input since specific format fields are 
not required. Variable names must be punched on the cards, 
thus allowing more rapid Identification of variables and their 
values. The use oi NAMELIST eliminates repetitious Input 

of data when running multiple configurations since only the 
variables with different values from the first case need to be 
Input. Note that this Is slightly different from the normal 
NAMELIST feature of updating the previous case. This Is done 
using temporary disk storage. Three NAMELIST blocks are 
used: 

MAINl includes all variables previously read in 
data groups In the MAIN program, except for NRUNS and the 
tabular Input of the propellant temperature distribution 
(ITEMP - 0). 

AREASl Includes all variables previously read In data 
groups In subroutine AREAS, except for the tabular input 
of burning areas. 

AREAS2 includes all the variables for the tabular input 
of burning areas (INPUT - 1 or INPUT - 3). 

21 
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Ah example of the use of these NAMELIST blocks In multiple con- 
figurations is given in Table III-l. For multiple cases, blocks 
MAINl and AREASl must be Included In all cases even If no variables 
are listed (See ABEASl block of second case in Table III-l) . 
NAMELIST Is a stanoard FORTRAN IV feature, and additional infor- 
mation can be obtained from any good reference in FORTRAN IV. 

3. The use of NAMELIST now allows the elimination of the input 
variable NTABY, which was necessary to insure the reading of all 
tabular input data cards before proceeding to the next data case. 
When data is read from NAMELIST, preceding unread data cards are 
skipped, thus eliminating the previous need to do this using 
NTABY. 

4. The input variable MU is eliminated since it is only used for 
calculating the initial Reynolds number, which is not used 
elsewhere. Obviously the calculation and output of initial 
Reynolds number are also eliminated. 

5. An output option (100) is included in NAMELIST MAINl to allow a 
reduction in output volume. 

The options for 100 are; 

0 for partial output (2 lines) 

1 for full output (4 lines) 

Tables III-2 and III-3 illustrate partial and full outputs, re- 
spectively, and include various other input options. 

6. An input variable (PATMI) for the initial atmospheric pressure 
(in psia) is included in NAMELIST MAINl to accommodate pressures 
other than standard (14,696 psia). 

7. Numerous repetitious computations are consolidated, and some 
program logic is optimized. This does much to reduce the 
required execution time of each configuration. 

8. The convergence criteria (GMAX and GMIN) on the chamber pressure 
iterations are slightly relaxed. This gives a reduction 

in execution time, while not adversely affecting the program 
results. 

9. In the output data, FAV and FVACAV have been replaced by AVRAD 
and AVEPS, which are the average nozzle throat radial erosion 
rate (in in/sec) and the average nozzle expansion ratio, re- 
spectively. 
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Table III-2. Sample of partial output (100=0) from the design analysis computer program. 
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Table III-3. Sample of full output (100=1) from the design analysis computer program, also showing 
Ignition transient, propellant tenq>erature distribution and inert weight results. 
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10. The special relationship used to calculate head end chamber 
pressure for combination circular-perforated and star grains 
Is eliminated. The same approximate relationship Is now 
used for all grains, because detailed examination of results 
over a period of several years shows the special relationship 
tends In general to underestimate Che head end pressure. 


These modifications have resulted in substantial savings In both the 
number of Input data cards (approximately halved) and execution time of 
the design analysis program. With a full output (100 ■ 1, Table III-3) 
the reduction In execution time Is approximately 20 percent, and 
more significantly, with a partial output (100 •> 0, Table III-2) the 
reduction In execution time Is approximately 40 percent, and the volume 
of output Is approximately halved. This reduction In execution time 
should, however, not be misinterpreted as a reduction In total computer 
time or compilation time. While the compilation time comprises, by far, 
the majority of the total computer time for a single run, and has not been 
reduced by this effort, the compilation time can be totally avoided by the 
use of an object deck, l.e., a copy of the program which has already 
been converted to machine language and requires no compilation time. Thus 
the reduction i.i execution time is most significant when using an object 
deck or when running multiple configurations. 

The entirety of these modifications have also been applied to the 
design analysis program with grain deformation (Section IV) along with 
additional modifications unique to the grain deformation analysis and the 
addition of a new erosive burning law options. These additional modifica- 
tions are discussed in Section IV. 

Monte Carlo Computer Program 

Since the Monte Carlo Program is composed, In part, of the design 
analysis program, most of the modifications described above can be 
utilized. Modifications number 1, 6, 7, 9, and lU are applied to the 
Monte Carlo program. Modification number 2 is also adopted insofar as the 
non-statlstical input variables are concerned. The modifications described 
In number 7 are extended to include that portion of the program unique to 
the Monte Carlo program. In addition to these modifications, the input 
option IRAND is eliminated so chat only the random number generator RAMDU 
Is used. The GAUSS (machine Independent) random number generator Is 
eliminated as an option since the true randomness of this generator as 
adapted into the Monte Carlo program is questionable. An example of the 
Monte Carlo input data is Illustrated In Table III-4. Note Chat a change 
In the statistical input data format is Incorporated, in conjunction 
with the NAMELIST modifications, to reduce the total number of input 
data cards by approximately one-half. An example of the Monte Carlo pro- 
gram output Is shown In Table III-5. 

As In the design analysis program, a considerable savings in execution 
time Is realized in the Monte Carlo program using these modifications. 



Table III-4. Saople Input for the Monte Carlo computer program 
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Table III-5. Sample output for the Monte Carlo computer program. 
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In order to assess the reduction In execution tlme> a sample set of 6C 
motors (30 SBM pairs) Is analyzed. The reduction In execution time Is 
approximately 20 percent. Unlike the design analysis program, the 
execution time of the Monte Carlo Program represents, by far, the majority 
of the total computer time, resulting In a considerable savings even 
though no significant reduction In compilation time has been achieved. 



IV. EFFECTS OF GRAIN DEFORMATION ON INTERNAL BALLISTICS 

Reference 2 explains how deformation of the propellant grain can 
affect the propellant burning rate and account for a significant portion of 
the "scale factor" between large motors and small ballistic test motors. 

The important factor is the tangential strain at the burning surface of the 
propellant grain. A method developed by Smith (Ref. 9) was used In Ref. 2 
to estimate the strains at various times for circular -perforated grains. 

The method was coupled with the design analysis computer program to analyze 
the effects on Internal ballistics. Comparisons were presented of Internal 
ballistic test data and theoretical evaluations with and without grain 
deformation for the Titan III C/D and Castor TX354-5 SRNs. Although good 
results were obtained with the grain deformation analysis, It was felt that 
additional verification of the method was needed before the analysis is 
Incorporated Into the Monte Carlo computer program. Also, the effects of 
star-type grain deformation needed to be Investigated. 

It was planned to test the hypothesis on circular-perforated grain 
deformation on two additional SRMs: the TU-A55.02 (64 in. dia.) and the 

Space Shuttle nozzle erosion test motor. However, the erosion test motor 
proved inappropriate because the grain was not case-bonded. Also, the 
anomalous performance of the TU-455.02 Indicated tliat additional factors 
were Involved which affected the burning rate (Ref. 10). The Irregular 
behavior of this SRM thus masked our early efforts to ascertain the effects 
of grain deformation. This led first to investigation of potential effects 
of nozzle throat variations during SRM operation and finally to a seml- 
cmpirlcal analysis of propellant erosive burning ch,aracteristlcs. Details 
of these analyses are presented In this Section along with comparisons of 
test data for the TU-455.02, Castor TX354-5 and the first Space Shuttle 
SRM (DM-1) with theoretical results using the grain deformation analysis 
and a new burning rate model. Also presented is an Investigation of the 
effects of star-type grain deformation. This Includes both an experimental 
test using an inert subscale model of the Space Shuttle star grain segment 
and a theoretical finite element analysis. 


Investigation of Nozzle Throat Variation 

Several large SRMs Including the previously mentioned TU-455.02 were 
shown by Ref. 10 to have actual pressure-time curves which did not match 
theoretical predictions. These motors were found to have a consistent 
pattern of deviation. The pattern was such that the actual pressure was 
lower than theoretical over the first and last portion of the trace and 
higher than theoretical during intermediate portions. The deviation was 
explained in terms of a burn-rate variation with web thickness. A tern 
known as the burn-rate anomaly rate factor (BARF) was defined as the ratio 
of actual burn-rate to theoretical burn- rate. Reference 10 discussed 
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some of the possible causes of the BARF variation but did not present any 
theory to completely explain the phenomenon. Some of the possible causes 
mentioned were grain design, deformation of grain, temperature gradient, 
and propellant characteristics. 

The design analysis program with grain and case deformation effects 
was used during the present study to verify Chat deformation could not by 
Itself account for BARF variations. Indeed, In the case of TU-45S.02, the 
deformation effect seems to Increase the Initial BARF. The negative result 
raises Che question of whether burn>rate is the parameter that varies with 
time or whether other parao»ters cause the burning rate anomaly. 

Notale throat variations during motor operation is next examined to 
determine if the deviation In chamber pressure could be caused by a noztle 
throat area deviation. The scheme used to deduce nosele throat diameter, 
in Che absence of direct measurement. Is to compare the experimental thrust 
and head end pressure for a particular SRM to obtain the apparent throat 
diameter. This gives throat diameter as a function of time. The noasle 
dlemeter can thus be found. Independent of any burn-rate effects. 

The first step In this Iterative procedure is to assume a noszle throat 
area. Next, the Mach number at the nozzle end of the port Is determined 
from the ratio of the area of the port to the area of the throat, using the 
isentropic relation. 


M - (A*/A){211-Ky-1)m2/2]/(y+1)}^'^'^^^^^^^''^^ . (IV-1) 

The nozzle end stagnation pressure Is found from the experimental hoad 
end pressure from the following equation which can be obtained from 
relationships given In Ref. 11, 

W^h " U-KY-l)M^/21^^^^"^V(l-HrM2) . (IV-2) 

Eq. lV-2 is based upon a constant port area. Applied to a quaal-steady 
state aicuaclon. It is Independent of the burning rate distribution over 
the propellant surface. 

The exit Hach number is determined from the nozzle expansion ratio 
using Eq. lV-1. The exit pressure is obtained from the exit Hach 
number and the stagnation pressure. Assuming isentropic flow in the nozzlo. 




(IV-3) 
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The ideal thrust coefficient is then computed from 
o 

The thrust Is obtained from 


0 

and Is compared with the experimental thrust F*. If the difference betyov>> 
computed and experimental thrust Is greater than a specified tolerance , a 
new A* Is computed from 


-ir 



A* - [F'/C„ - A (P - P )]/(P Cp ) (lV-6) 

F e e a on F ' 

o 

and the Iteration continues with a new evaluation of M ; otherwise the tiir. 
and port radius Is Incremented and the next experimental values used. 

The SRMs Investigated using this procedure are the TU455.02, the 
LPC-156 Jet Tab, the AGC-260-1, the AGC-260-2, the TCC 156-2C-1 and the 
TX35A-5 of References 12 through 17, respectively. 

The variation In nozzle throat diameter with time for each of thest 
motors Is given In Fig. IV-1. A linear Interpolation between the actual 
Initial and final throat diameter Is also shown on the figure. Some error 
In the results can be attributed to the necessity of obtaining data from 
plots In test reports. Also, for several of the SRMs which were fired 
vertically. It was necessary in the test reports to add the total weight «>i 
propellant burned at each time to the corresponding thrust. Inaccuracy in 
estimating the burned mass can also contribute some error to the results. 
However, the motors generally show an Initial apparent closing down of tho 
throat possibly due to mass blocking effects from low available energy 
propellant gases emanating from the. aft end of the grain. The analysis 
does not provide a scheme to account for such two-dimensional nozzle flov 
phenomena. An additional source of error in the analysis is that some pt 
the SRMs have a slight axial taper to their ports while Eq. IV-2 is 
based on constant port area. However, burning rate is used in the analysiN 
only for evaluation of the port diameter as a function of time, and errors 
from this source should be small at least for the initial portion of the 
trace where there appears to be some effect of nozzle throat variation not 
accounted for In the usual nozzle erosion models. Beyond this, the results 
are Inconclusive and do not by themselves explain the burning rate atwnui ly . 




Figure IV-1, Variation In 
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Burning Race Model 


Since the nozzle throat area investigation failed to give a clear 
explanation of the burning rate characteristics of the SRMs examined, 
attention was directed to other possible causes of the anomalous behuvloi . 
The TU-A55.02 was selected for the Initial investigation because the BARK 
effect was so pronounced in this SRM. Although Ref. 10 examines a variety 
of potential causes of the burning rate anomaly, no investigation of the 
contribution of erosive burning is discussed ■'n the reference. This Is 
possibly because measured chamber pressures are generally lower than pre- 
dicted for the initial portions of the pressure-time traces in- the TU-A53.02 
and other SRMs examined In Ref. 10, and this is the reverse of what one 
would expect from erosive burning. However, we observed that the measured 
chamber pressures are initially low in the. motors exhibiting BARF and it 1:- 
well known that erosive burning rates are quite sensitive to pressure 
(Ref. 18). This led to an investigation of erosive burning as a possibl.- 
cause of the burning rate anomaly. 

The first effort was strictly empirical and was directed toward 
determining a coefficient, a(u), and an exponent, n(u), in the relationship 


r 


a(u) r 


,n (u) 


(IV- 7) 


that gives a good match between test data and theoretical performance. The 
design analysis computer program was used to obtain the predicted perfonwanco. 

The design analysis program calculates two burning rates: one neai 

the nozzle end of the propellant grain at the point of imixlmum Mach numbti 
and one at the head end where u“0. To obtain the total mass flow rate 
upstream of the point of maxltium Mach nvimbor a linear average of the t\%’o 
burning rates is taken (Ref. 5). Wlille introducing some question as to the 
accuracy of the resultant calculations, the use of this simplified progi am 
permits economical assessment of SRM perfornuince (Over 100 sopaiate pei 
formance predictions were made in the Investigations discussed in this suh- 
section). The appropriateness of linear averaging of burning rates is 
discussed furt'.*er later in the report and a modification to the method intro- 
duced. In calculating the head end burning rate, the a and n corresponding 
to small ballistic test motor rates, assumed to give zero velocity behavior, 
were used wherever possible for each SRM examined in the course of this .stndv 
so as not to confound the investigation with artificial "scale factors." Kin 
the TU-455.02, it was i\ccessary to estimate, the zero velocity rate from 
performance of the large motor somewhat prior to talloff since the ballisiii 
test motor rate was unavailable. Grain deformation was not accounted for in 
the first evaluation. 

After a large number of runs with the design analysis program wherein 
various simple functions of u for a and n at the nozzle end of the grain 



35 


were tested, a fair match was obtained. The results were not Improved In 
general by consideration of grain deformation. Similar empirical express 
slons failed badly when tested for Castor TX354-5. The Investigation did 
reveal that there is an apparent general depression of the burning rate 
coefficient and increase in burning rate exponent in Eq. IV-7 at the 
higher flow velocities which is not inconsistent with the qualitative erosive 
burning behavior of many propellants. In partlrului.' it has been demonstrated 
(Ref. 18) that some propellants exhibit a threshold velocity below which the 
burning rate is lower than for zero velocity. The threshold is In general 
lower for higher pressures and one might alternatively speak of a threshold 
pressure below which velocity depresses burning rate. The threshold velocity 
(or pressure) is not usually very high but it should be kept in mind that 
even for an SRM with small port-to-throat area ratio a substantial portion 
of the propellant surface is subjected to low velocities even during the 
initial portion of SRM operation. Thus it seems that the low velocity 
regime is worthy of additional investigation. Several of the more widely 
used erosive burning rate laws do not even permit representation of burning 
rate depression, notable among which is the law of Robillard and Lenoir 
(Ref. 18). 

Since a relationship of the form of Eq. IV-7 having general applicabillLy 
was not found, the investigation turned to an examination of alternative 
relationships. Poor results have been oh alned in previous investigations 
with the design analysis program using tae law of Robillard and Lenoir 
(Ref. 5). This law uses an approach which is more or less conventional in 
its treatment of the erosive burning effect as one that is strictly additive 
to the rate obtained at zero velocity. In the present investigation, we take 
the alternative approach of considering the erosive effect to be one of 
direct reduction of the flame height and, following Summerfleld (Ref. 18, 
p. 38) on ammonium perchlorate composite propellants, take the burning rate 
to be inversely proportional to the average flame height. 

The usual form of the Summerfleld formula applicable in the absence 
of erosive burning is 


1/r^ - c^/P + Cj/P^/^ 


(IV-8) 


The formula fits experimental burning rate data about as well as ro ■ ap 
(Ref, 18, p, 38). The first term on the right-hand side of Eq. IV-8 
represents the contribution of the diffusion flamelets to the average flame 
height while the second term gives the effect of the surface covering by 
premixed flames. In the presence of erosive burning, Barrfire and Larue and 
others treat the constant c, as a function of velocity and the exponent of P 
in the last term of Eq. IV-8 is assigned a value of 2/3 (Ref. 18, p. 436 
and p. 445). Our results with attempts to match test perfornumce of Tl'- 
455.02 with similar representations are unsatisfactory. 
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We next hypothesize that a reduction In flame height and a correspond- 
ing Inversely proportional Increase in burning rate accompanies the 
establishment of a flow velocity over the surface of the propellant and 
the effect on flame height Is subtractive so that 


1/r = 1/r^ - C3ReO-5(i_p^^/p)/L (IV-9) 

Here r^ =» aP” Is used to represent the burning rate in the absence of 
velocity rather than the more complicated relationship, Eq. IV-8, as a 
matter of convenience. The form for the last term in Eq. lV-9 evolved 
after a number of comparisons with the test performance of both TU-455.02 
and Castor TX354-5. In Eq. IV-9, Re is the flew Reynolds number and it was 
found during the comparisons that best results are obtained with a Reynolds 
number based on flow-path length. Also, it was observed that the ratio of 
the pressure above a certain critical pressure P^r to the pressure P itself 
plays an apparently important role in the reduction of flame height over 
and above that played by Re, hence the term (1-P(,j./P). The investigation 
was conducted to obtain the best matches between theoretical and test per- 
formance without grain deformation being taken into account. About equal 
quality results are also obtained for the two SRMs with an exponent on Re 
of 0.8. However, when grain deformation effects are included, the exponent 
of 0.5 proves better. Higher powers of (1-Pj.y/P) were also considered but 
give poor results. 

The viscosities of the propellants in this st^.dy are considered to be 
substantially the same so that the final form of the equation used to 
calculate the burning rate at the point of maximum flow Mach number is: 


1/r = 1/r^ - c^(GL^^^/L)0-S (1-P^^/P) (JV-30) 

Here L^ef is the initial "controlling" grain length of the TU-455.02 
(205.75 in.) which is used as a reference dimension in the development and 
G is the mass flow rate per unit area at tlie point of maximum Mach number. 

Although results with the modified flame height law (Eq. IV— 10) com- 
bined with the grain deformation are quite good for the TU-455.02 and 
TX354-5, for both SRMs, the theoretical analysis yields pressure- time 
traces that rise progressively, somewhat above the actual performance traces, 
as tailoff is approached. One explanation is that Eq. IV-10, while con- 
taining a threshold pressure, does not represent the velocity sensitivity v>f 
this threshold; that is, the threshold pressure la the same at all velocities 
so that as long as P > P^r* the erosive rate is greater than the zero 
velocity rate no matter how low the flow velocity. This discrepancy is 
adjusted to some extent in the new analysis by modifying the distribution of 
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burning rate over the grain. In place of using a linear average of head 
and nozzle end rates to calculate the mass generated, the average burning 
ratio is taken as 

r^^ - 0.5[l+(y/b)5]r^ + 0.5 [l-(y/b)5]r^ (IV-11) 

where b represents the average web thickness of the main propellant grain: 

b - - |z/2l - (y<b) 

b - y (y>b) (IV-12) 

These relationships cause the average burning rate to rapidly approach the 
zero velocity rate in practical SKMs as talloff is approached while having 
negligible Influence during about the first two-thirds of web action time. 
Although strictly empirical, the approach seems reasonable for at least 
partially compensating for the fact that the linear average clearly becomes 
Inadequate as burning progresses and more and more of the burning surface 
is subjected to less than threshold velocities. 

Final Comparisons with Test Data 

The final comparisons with test data are presented .'<n Figs. lV-2 and 
IV-3 for the TU”455.02 and TX354-5, respectively. Results are given 
based on the modified flame height model both with and without grain de- 
formation being taken into account. Values of the constant material 
parameters used are identified on the figures. It is notable that the best 
fits are obtained for both SRMs when the grain deformation effect is in- 
cluded in the evaluations. 

As a final test, an evaluation was performed for the Space Shuttle SRM 
and comparison made with the test data from the first development motor 
(DM-1) using the "as built" configuration. The Space Shuttle SRM has a 
FBAM type propellant quite similar to that of the TU-455. 02 except that a dif- 
ferent curing agent and a somewhat different percentage of burning rate 
catalyst is used. For the Space Shuttle evaluation, the TU-455. 02 constants 
were used in the flame height law. For the Space Shuttle deformation 
analysis, the reported Poisson's ratio of 0.499 was used while for the 
TU-455. 02, the reported ratio of 0.4995 was used. The results are shown 
in Fig. lV-4. In this case, a better match is obtained without the grain 
deformation. However, the analysis with grain deformation would be as good 
or better if a Poisson's ratio of 0.4995 were us> d in the Space Shuttle 
evaluation and would be definitely better than the corresponding analysis 
without deformation if the erosion constants were adjusted for the best fit 
with grain deformation using the reported values of Poisson's ratio. 

We feel that the results of this study provide additional support 
for the validity of the '^rain deformation analysis as well as a new approach 




figure IV-2, Comparison of test data for the TU-455.02 SRM with theoretical 
perfonnance calculated using the modified flame height burning 
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Figure IV-4 . Comparison of cest data for th Spa;e Shuttle SRM (DM-1) with thee 
performance calculated using the modified flame height burning rat 


to representation of erosive burning phenomena that is worthy of further 
examinatior. However, the results for the Space Shuttle SRM introduced 
some questions about the applicability of the grain deformation that 
should be resolved before this approach is incorporated into the Monte 
Carlo program. Also, the burning rate model should be tested in a more 
comprehensive analysis which computes local burning rates at a large 
number of stations throughout the propellant length while providing for 
threshold velocities in a more realistic manner than the simplified design 
analysis program. Improvement of the basic burning rate law expression 
is clearly possible. For example, a brief investigation has shown that 
the law as presently formulated does not fit the performance of small 
ballistic test motors but tliat betror results are probably obtainable with 
larger exponents on the Reynolds number term. 

It also may be possible to find a better approach than tiat provided 
by Eqs. IV-11 and IV-12 to account for threshold velocity in the simpli- 
fied design analysis program. The fact that, for each of the three SRMs 
examined, the predicted traces still tend to rise above those of the 
measured data as tailoff is approached is quite possibly due to inade- 
quacy of the design analysis program in its artificial treatment of the 
threshold velocity effects. 

Further study of the effects of nozzle throat area variations may 
provide a means of improving prediction or chamber pressures for the 
initial portion of the trace where, fo’' the three cases presented, the 
predicted pressures are low. This phenomenon as previously noted might 
well be due to mass blocking of the ;hroat by low available energy flow 
from the aft end of the propellant. Simple representation of this 
blockage may be possible within the framework of the design analysis 
program, and, of course, would be applicable to more sophisticated pro- 
grams employing the modified flame height model. 

Deformation of Star C ins 

A portion of the Space Shuttle SRM head end segment contains a star- 
type grain. In the grain deformation analysis using the modified design 
analysis program, deformation of star grains is disregarded. We next 
examine the validity of this assumption making use of both experiment 
and theory with particular emphasis on the Space Shuttle SRM. 

For the experimental Investigation, a scaled model of the cross-section 
of the forward segment of the Space Shuttle's Solid Rocket Motor was con-* 
structed. An aluminum mold was constructed with dimensions that were 1/J6 
of those of the actual Shuttle motor. 

The material selected to represent the grain was Dow Coming's 3110 RTV 
silicone rubber. This material consists of a white liquid encapsulant and a 
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curing agent. The cured material has mechanical p:*operties approximately 
equivalent to those of the Space Shuttle propellant. Uslr'\ a mii.tng ratio 
of 1 part curing agent to 10 parts of encapsulant, the lu.txturo was prepared 
and poured Into the aluminum mold. The silicone rubber was allowed to cmo 
for 48 hours. The model was then separated from the mold and the excess 
material was trimmed off. The thickness (length) of the finished model 
1/2 In. and the diameter 4 In. 

A three-piece chamber was constructed from a 1/2-inch thick sheet of 
Plexiglas (See Fig, IV-5), The center portion of this chamber consisted 
of a Plexiglas ring to which the silicone model ' ■ glued to represent case- 
bonding. The ring was then sandwiched b< tween t’ Jexlglas discr ana the 
three pieces were bolted together with Ixteen bOj.cS. Clearance space was 
provided on both sides of the model and the Plexiglas discs. A pressure tap 
was Inserted through the center of one d.isc. A steel cylindrical contain- - 
filled with mineral oil was connected by hoses to the Plexiglas chamber 
A steel plunger was Inserted into the top of the cylindrical container so 
that a load could be applied to the plunger and model by a universal tcru.lon 
compression testing machine. A Bourdon-tube pressure gage was connected to 
the hose between the cylinder and the piston (top photograph, Fig. TV-5), An 
extensometer was mounted as shown in the bottom photograph on Fig. iV-S f< r 
the purpose of recording changes in diameter of the chamber. 

A five-inch rectangular bar was poured from the same mixture tiiat wa.-' 
used to cast the model. After curing, this bar was inserted into the univi-isa 
testing machine. At various loads, the length and width of the bar were 
measured. From a plot of the stress-strain curve determined from the elniigii 
tlon of the bar as a function of applied load, the value of Young's 
was found to be 500 psl. By plotting the change In width against the ciuu ge 
In length at each load, the value of Poisson's ratio was determined fiom ( he 
slope of the resulting curve to be 0,498. 

A camera with a close-up lens was pieced in front of the assembled el.nihe 
and model. The camera was focused so that one star point com’'' etely lilled 
the frame. A picture was taken at 0 psl (gage pressure), Us.ng the nulv. -.i'. .il 
testing machine, a load was slowly applied. A picture was taken at each 
psl increment up to and including 400 psl. The photographs were developed 
as 8 X 10 enlargements. By tracing the perimeter of the star point at 0 psl 
and 400 psl, the change in perimeter could be calculated. xYlso, by super- 
imposing these two traces as in Fig. IV-S, the deflection of tlie stai p>'!i>: 
could be observed. Figure IV-6 shows the deflection of the star point .it 
400 psi as compared to its initial position at 0 psi. The point was totmd to 
elongate slightly in the radial direction and contract in the c Iroumteiont lal 
direction. The perimeter was found to Increase by 0 , 22 % at 400 psi. 

A theoretical Investigation of the dofcrmatlon ef star grains Juo to 
pressurization was also performed using a finite element comput<‘r pi >gi.".. 

The purpose of this investigation was to determine the magnitude of defi'ini.i- 
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tion, correlate the results with experimental data and, If necessary, apply 
proper adjustments to the design analysis program, preferrably with a simple 
approach. 

The computer program used In this analysis was obtained from Ref, 20 
and was modified for use on the IRM 370/1S8. The program performs a two- 
dimensional stress-strain analysis of a propellant grain, for a condition of 
either plane strain or plane stress. A condition of plane strain was used 
for this Investigation since as shown In Ref. 9 when the grain Is case- 
bonded longitudinal strains arc small. 

The accuracy of the program was tested using a finite element model of 
a circular-perforated grain. Reaulti; of this test Indicated that deformationa 
computed by the program were consistent with those computed with the equations 
of Refs. 2 and 9. A major subprogram was added to the finite element progiam 
for automatic grid generation of finite element models of standard and 
truncated star grains using inputs identical to the design analysis program. 
Another subprogram was coded for computing the change In perimeter of standard 
and truncated star grains due to deformation. 

The Space Shuttle SRM star grain initial geometry was Input to the 
finite element program with a cluimber pressure or 900 pslg. The computed 
change in perimeter was -5.3 In. or a decrease of 0.46J. This amounts to a 
total change in burning area of -939.6 in.^ for the star grain or a decrease 
of 0.212 in total burning area including the circular -perforated grain area. 

A finite element model corresponding to the subscale experimental model was 
also Input to the program with a pressure of 400 pslg. The cose modulus 
value was adjusted to obtain theoretical case deformations eqvilvnlent to the 
measured case deformation. This resulted in a decre.ase In perimeter of 0.3* 
as opposed to the increa.so of 0.222 measured In the experiment. 

For this particular grain conf Igur.at ion the controlling pamraeter tor 
the sign of the perimeter change is the length of the st.ir point as compared 
with the web thickness and the width of the star point. Although the experi- 
ment was designed to approximate a condition of plane strain on the model, 
end effects may explain the differences. For the experimental model, 400 
pal pressure was presunuibly applied to the end faces of the nwdol, whorcas, 
the finite clement plane strain solution vlelded axial stresses on the 
ends which were 0 to 40 psi less depending »... the radial po.slt ion. The 
lower end pressures would clearly cause the sla. -'olnts to he shorter and 
could account for the dneroase In perimeter In the thcoi'ct leal !’.ol\Ulon. Also, 
the finite clement soletion mokes use of elements with straight sides pro- 
ducing some inaccuracies which cun bo minimised by Increasing the numlnn' of 
elements. However, for the pavtlrvilnv configuration tested, regardless of 
the sign, the magnitude of area change is snwtll with roapect to that pvodviced 
by the circular-perforated grain dofonmttlim. 

For other configurations of greater web thicknesses and shorter star 
points, the effects nwy be more significant . In general, the finite element 
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program could be used to relate pressure and star grain configuration 
parameters to change In perimeter. However, an overall program for Internal 
ballistic analysis would be somewhat cumbersome since the perimeter changes 
would have to be evaltiated at each Increment of time. He feel that at least 
for star grains similar to those of the Space Shuttle, the effect of star 
grain deformation is Insufficient to warrant the application of such a 
program. The studies thus far completed could prove of considerable value In 
finding a simplified approach to analysis of star grain deformation effects 
that would be consistent with the economy objective of the design analysis 
program. 


Modification of the Deformation-Design. Analysis Computer Program 

All modifications to the design analysis program of Section III have 
been Incorporated Into the grain deformation design analysis program of 
Ref. 2. The purpose of these modifications as mentioned In Section 111 Is 
to reduce computer execution time, program size, and Input and output com- 
plexity. The equations for the modified burning rate model discussed In 
the present section have also been Incorporated Into the design analysis 
program as an option. 

The basic Inputs to the design analysis program with c.p. grain de- 
formation and modified burning rate model options are Identical to those 
for the design analysis program of Section 111. Also required as Inputs 
tinlque to the new program are the following variables specified In NAMELIST 
MAINl: 


1. Deformation Inputs (Do not use If GRAIN ■■ 2) 

ISO ■ 0 for no deformation effects 
■ 1 for deformation effects 
PItt)D - elastic modulus of the propellant, psl 
CMOD ■ elastic modulus of the motor case, psl 
FMU <■ Poisson's ratio for the propellant 
CMU ■■ Poisson's ratio for the motor case 

ALPTS ■ Linear coefficient of thermal expansion of the propellant, 
In/ In 

TAUC - Motor case thickness, Inches 
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2* Hodlfled burning rate model inputs 

IHO •* 0 for al^ burning rate law 

• 1 for Roblllard-Lenoir erosive burning law 

■ 2 for modified flame height model (Do not use if IGO ■ 1) 

ERCON ■ Coefficient of the erosive burning term for the modified 
flame height model, (sec^ /slug) 


PCRIT ■ Critical pressure in the modified flame height model, psla 
LREF ■ Reference length In the modified flame height model. Inches 

It Is not necessary to Include these new variables when the effects 
are not desired; that is, when IRO or ISO Is zero. Sample Inputs and 
outputs for the grain deformation design analysis program with the modified 
burning rate law options are Illustrated in Table IV-1 and IV-2, respectively. 
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V. THERMOMECHANICAL ANALYSIS 


In Ref. 2. the variation in propellant burning rate due to thermo- 
elastic coupling was shown to be a function of both elastic modulus and 
pressurization rate. Effects of the thermoelastic coupling on propellant 
burning rate was shown to be small in Ref. 2. However, the propellant 
modulus was chosen constant and as such was independent of spatial varia- 
tions in temperature and strain rate. Both temperature and strain rate 
are known to have significant influence on the propellant modulus. Hence, 
it was decided to improve the propellant modulus model to account for 
spatial variations in temperature and strain rate before making a decision 
on the Influence of the thermoelastic coupling. 

In Ref. 2 an exponential pressure rise was used for the mechanical 
load. Under certain operating conditions the propellant may be subjected 
to an oscillating pressure load. The Inherent irreversibility of the 
thermal processes produced by thermoelastic coupling makes this cyclic 
loading condition a likely candidate for producing changes in propellant 
temperature and hence in propellant burning rate. This loading condition 
is investigated in the present work using the improved model for the 
propellant modulus. The effects on thermoelastic coupling of varying 
Poisson's ratio, thermal conductivity and coefficient of expansion are also 
Investigated in the present work. 

A brief analysis of the effects of mechanical vibration on motor 
performance, including energy dissipation associated with viscoelastic 
materials is presented. Also discussed are the results of an Internal 
ballistic analysis of ignition transients incorporating the Zeldovlch 
and Novozhilov (Z-N) transient burning rate model (Refs. 2 and 7). 

Propellant Modulus Model 

The strain rate and temperature dependency of the modulus were modeled 
using the relaxation modulus model of Blatz (Ref. 21) shown in Eq. V-1. 


and the well known time- temperature shift factor of Williams, Landel and 
Ferry (Ref. 22) given by 


log 


10 T 


c;,(T. 


■T^)/(f:2+T-T^) 


(V-2) 
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where the time, t, is given by t s.- At/a^ and At is the time step used in 
the numerical solution of the thermoelastic problem. 


The constants 3^ end in Eq. V-1 and the constants C£ end 

Tg, the glass transition temperature in Eq. V-2» were obtained from 
Ref. 23 for Thlokol propellant TP-H8156. The constants Eg, the steady 
state or long time modulus, and Eg, the modulus at the temperature Tg, 
were adjusted to obtain a value of approximately 550 psl at 70°F. 

This was done because the constant propellant modulus used in the base- 
line analysis of Ref. 2 was 550 psl at 70*F. 

The values used for the constants in Eqs. V-1 and V-2 are given in 
Table V-1. 


Table V-1 


Propellant Modulus Constants 


3 0.164 

o 


2.89 X 10~® sec. 


E 100 psi 


31,25C p^l 
-17.44 ^ 


51.6*K 


243.15*K 


A plot of the modulus versus temperature is shown in Fig. V-1. 


The model described above was Incorporated into the existing thermo- 
elastic analysis computer program. An exponentially increasing pressure 
load, identical to that used in Ref. 2 was used with the new propellant 
modulus model. The results obtained are shown in Fig. V-2 compared to 
the baseline model in Ref. 2. In Fig. V-2, the ratio of the dynamic burn- 
ing rate, r, obtained using the relaxation modulus and the Z-N model of 


Refs. 2 and 7 to the "steady state" burning rate, r , 

85 


obtained from 


(V-3) 



PROPELLANT 



Figure V*l. Relaxation modulus versus temperature. 
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Is plotted versus time. The results using the new modulus model show 
essentially no difference between the two models except during the Initial 
portion of the curve where there Is less Influence of thermoelastic 
coupling than was obtained previously. This Is due to the lower value 
of propellant modulus at the higher temperatures near the burning surface. 
These results Indicate that effects of thermoelastic coupling are probably 
negligible under the conditions anticipated In a Space Shuttle type SBM. 

Figure V-3 shows the effect on propellant burning rate caused by 
changing the Initial pressurization rate from 14000 to 28000 psl/sec. The 
relaxation modulus was used for both cases shown and the results indicate 
the same trend established In Ref. 2 of a higher pressurization rate pro- 
ducing a larger effect. These results are presented for the condition which 
Includes thermoelastic coupling. However, the omission of the thermo- 
elastic coupling produced no significant differences. The results show 
that the dynamic burning race approaches the steady state burning rate very 
ri^'idly and at 50 milliseconds there Is less than 4 percent difference for 
b- th pressurization rates. The pressure- time history for both loading 
conditions Is also shown In Fig. V-3. 

Influence on Propellant Burn Rate Produced by Oscillating Pressure Loads 

The propellant burning rate Is directly proportional to the pressure 
In the combustion chamber; hence, If the combustion chamber pressure 
oscillates, the propellant burn rate also oscillates. Perturbation and/or 
oscillation of the combustion chamber pressure may be produced from many 
sources. A few sources of particular Interest are wave motion In the 
combustion chamber during Ignition or rap . thrust termination, irregular 
and possibly unstable burning, motor and/or vehicle vibrations and in the 
case of vehicles like the Space Shuttle, vibrations produced by unsteady 
combustion in the liquid engines; e.g. , the POGO phenomenon. The latter 
problem Is studied In Ref. 24 and was shown to be Insignificant with regard 
to the Space Shuttle propellant. The effect of motor and vehicle vibrations 
Is discussed In a later portion of this section. For the present we con- 
sider osclllatlcns In combustion pressure such as might be produced by a 
aeries of equal amplitude pressure waves moving across the surface. The 
pressure oscillations are represented by functions of the following form: 


P « c, P sin (lit + P (V-A) 

1 av av 

where P^^ is the average pressure and Cj^ and co are arbitrary constants. 

This pressure function was Incorporated Into the thermoelastic 
computer program as the driving mechanical load with P^^, la and being 
used as variables. Samples of the results obtained are shown in Figs. 

V-4 through V-6. 
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Figure V*4, Results obtained for oscillating pressure loads 
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Figures V>4 and V-5 show the results for values of ranging from 
200 to 500 pal with u > 200ir rad/sec and Ci • 0.10. These results again 
show that only a negligible effect Is produced by thermoelastic coupling 
and hence only one curve, that which Includes thermoelastic coupling, Is 
shown. The ratio r/rgg Is plotted versus time and the curves Indicate a 
general tendency for this ratio to oscillate at a relatively small con- 
stant amplitude at low pressures but at higher pressures there appears 
to be a tendency for the oscillations to Initially be relatively large 
and then to decay. The oscillations have a period equal to the period 
of the pressure oscillations but are shifted slightly In phase. Possible 
causes for the phase shift may be the finite time required for the pro- 
pellant to respond to the pressure loading or the finite time step used 
In the numerical calculations; however, neither explanation has been 
verified. 

Figure V-6 shows the results obtained when the frequency for the load- 
ing condition of Fig. V-5b Is doubled. This was done to amplify the thermo- 
elastic coupling. However, the results indicate the thermoelastic coupling 
Is still negligible and the tendency for the oscillations to damp is still 
present. The reason(s) for the obvious change In frequency of the curve 
shown In Fig. V-6 are not presently understood. 

In summary, the oscillating pressure load showed no tendency to en- 
hance the effect of thermoelastic coupling phanomena. Therefore, It must 
be concluded that even though thermoelastic coupling does exist and can 
be modeled at least in an approximate way, there appears to be no advantage 
in considering Its effect on the burning rate of solid propellants under 
the loading conditions analyzed In this work and in Ref. 2. 

Effects of Poisson’s Ratio, Thermal Expansion and Thermal Conductivity 

In addition to the propellant modulus the energy equation for the 
coupled thermoelastic problem, 

~ - r(t)3T/3x + (A/;c) v2t - ToEe/[pc(l-2v) ] (V-5) 


contains terms involving other material properties. Three of these are 
Poisson's ratio, v, coefficient of thermal expansion, a, and coefficient 
of thermal conductivity, A. The propellant density, p, and specific heat 
at constant volume, c, also appear but are subject to less variation than 
the other three. 

The thermoelastic dissipation term, D, 


D - TaEe/[pc(l-2v)] 


(V-6) 



varies Inversely with Poisson's ratio. For a typical solid propellant, 
Poisson's ratio Is of the order of 0.499. With this value the term 
(l-2v) In the denominator Is 0.002. For a Poisson's ratio of say 0.4, 
the term (l-*2v) becomes 0.2. Therefore, a change of approximately 
25 percent In Poisson's ratlc changes a factor In the denominator by two 
orders of magnitude. Hence, the already small dissipation term becomes 
negligible with the smaller Poisson's ratio. It Is necessary tc conclude 
that reducing Poisson's ratio with all other factors unchanged essentially 
eliminates the thermoelastic coupling. By the same token. Increasing 
Poisson's ratio to say .4999 will Increase the coupling by an order of 
magnitude. Figure V-7 shows the results obtained for values of Poisson's 
ratio of .499 and .4999. Values for smaller Poisson's ratio approximate 
closely the uncoupled solution which has already been shown to be un~ 
Important. Again, only small changes result when the value of Poisson's 
ratio is changed. 

The remaining two material properties, coefficients of thermal con- 
ductivity and thermal expansion were varied Individually. The results 
obtained from this parametric analysis are shown In Figs. V-8 and V-9. 

Figure V-8 shows the resulting r/fgs values of the coefficient 

of thermal expansion of 0.1 x 10"^ and 0.1 x 10~‘'/"R. Only small changes 
In amplitude exist between the two cases with no apparent phase shift In 
the time-dependent response. Note also that the amplitude of t/r^g is 
essentially the same as that shown in Fig. V-5b using the nominal value for 
the coefficient of thermal expansion. Hence, it appears that variations 
In the coefficient of thermal expansion will produce only small changes In 
the dynamic burning rate. It should also be pointed out that the results 
shown In Fig. V-8 Include thermoelastic coupling, t^ut no significant effects 
were noted due to the coupling. 


Figure V-9 shows r/rg 3 values for coefficients 
of 5.1 X 10-5 and 5.1 X lO"'^ BTU/sec-ln.-*R. 


. thermal conductivity 


In these results there Is a significant change In amplitude as well 
as a phase shift produced by varying the thermal conductivity The lower 
thermal conductivity produces smaller amplitudes of r/rgg which means that 
the dynamic burning rate Is closer to the steady state rate than the dy- 
namic burning rate for the higher coefficlenc of thermal conductivity. 

T*’ere Is also a phase shift between the two curves due to the lower thermal 
conductivity slowing down the heat transfer process. This Is evident since 
the pressure loadings In each case were In phase, yet the peak amplitudes 
for the case of the lower coefficient of thermal conductivity occur at later 
times. Again, no significant Influence of thermoelastic coupling was noted. 

In summary, the results obtained from varying the propellant properties 
show that the propellant burning rate has a small variation with respect to 
these material properties. The burning rate variations produced by 
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Figure V-9. Comparison of results obtained for two different coef- 
ficients of thermal conductivity. 

(Pressure ■ 50 sin (200nt) + 500) 
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changing Poisson's ratio and the coefficient thermal expansion are In- 
dicative of modifications made to the thermoelastic dissipation term of 
Eq. V-5, since they only appear In this term. This Is also true for 
changes In the propellant modulus discussed earlier. The burning rate 
variations due to changes In thermal conductivity do not directly Indicate 
the Influence of thermoelastic coupling since the conductivity does not 
appear In the dissipation term of Eq. V-5. 

Effects of Mechanical Vibration on Motor Performance 


Mechanical vibrations In the propellant of an SRM produce temperature 
changes In the propellant due to viscous dissipation of mechanical energy 
as well as from volumetric changes associated with the theruioelastlc 
effect. The previous subsection has shown that the thermoelastic effect 
is negligible and it will be neglected In this portion of the work. 

A vibration analysis of a Space Shuttle type SRM is presented in 
Ref. 25. The analysis in Ref. 25 uses the results of Ref. 26 to determine 
the temperature change in the propellant per cycle of osci iatlon. 

For an element of propellant undergoing a pure shear vibration, the 
temperature change per cycle of oscillation determined in Ref. 26 Is 


AT = ttS jx 1 ^/G'pc 


(V-7) 


where AT is the temperature change per cycle, 3 is the shear loss 
modulus, G' is the real part of the complex shear modulus at a given 
frequency of oscillation, p is the propellant density, c is the specific 
heat at constant volume and x is the peak value of the oscillatory shear 
stress. 

Reference 26 also gives the temperature change per cycle of oscilla- 
tion for a uniaxial stress oscillation as 


AT - 3itBlol2/4E'pc 


(V-8) 


where 3 is the tensile loss modulus, E' Is the real part of the complex 
tensile modulus and a is the peak value of the oscillatory normal stress. 
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For a typical propellant we have from Ref. 24: 
6 - I - 0.4 


E' « 3000 psi 
G* ■ 1000 psi 
pc - 160 lbf/in2/*R 


These values give: 




AT 


uniaxial 


2 X 10-^ jap 


Since the temperature change due to shear oscillations Is larger, we 
conservatively assume that It holds for any vibrational mode. Reference 
25 gives 10 psi as a moderate to large oscillatory stress magnitude, using 
this stress magnitude, Eq. V-7 predicts that the change In temperature per 
cycle of oscillation Is 0.0008*R. Reference 25 also points out that exper- 
imental results from a vibration test of a Sparrow motor at Allegany Ballis- 
tics Laboratory Indicated a 10°R change after 30 minutes of vibration 
testing. The amplitude and frequency of the vibration test were not 
readily available and obviously the Space Shuttle SRM and che Sparrow motor 
are not of the same size class. But, If one makes the assumption that the 
Sparrow tests are Indicative of the magnitude of the temperature changes 
one could expect from a Space Shuttle SRM the value of 0.0008°R/ cycle would 
be considered a conservative estimate. With these limitations an analysis 
was made to determine the influence on performance of a Space Shuttle type 
SRM subjected to a vibrational environment which produced a 0.0008“R/ cycle 
change in grain temperature. The motor was assumed to vibrate with a stress 
amplitude of 10 psi at frequencies of 15 and 75 Hz. These frequencies are 
representativr of the fundamental vibrational frequency of the vehicle and 
of Its 5th harmonic. The results from this analysis for motor thrust 
versus time are shown In Fig. V-10 compared to the case of no vibration. 
Virtually no difference in performance Is noted between the no-vlbration 
analysis and the 15 Hz analysis. A small difference between the no- 
vibration and the 75 Hz analyses Is shown to exist. It Is Interesting to 
note that these differences occur Immediately preceding star grain burn out 
(- 20 sec), during the time of maximum dynamic pressure (- 60 sec), and 
during tailoff. By itself the vibration problem poses no real problem 
with regard to motor performance. However, due to the occurrence at 
these three critical times during the motor's operation, there exists a 
possibility, however small, of the enhancement by vibration of other per- 
formance variations to the point at which the combined effect could be 
important. The percent differences produced by mechanical vibration In web 
time were 0.11 and 0.56 percent for the vibrational frequencies of 15 and 
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75 Hz, respectively. The change in burning rate at web time was 0.19 
and 0.98 percent for the 15 and 75 Hz frequencies, respectively. The 
maximum head end pressure changed 0.05 and 0.23 percent and the average 
nozzle stagnation pressure changed 0.14 and 0.74 percent for the 15 
and 75 Hz frequencies, respectively. 

The results of this analysis indicate that the mechanical vibrations 
of the amplitude and frequency, 10 psi and 15 Hz, expected for a Space 
Shuttle type SRM would produce negligibly small effects on performance. 
However, it should be noted that larger amplitudes and/or frequencies of 
vibration tend to produce more significant effects on SRM performance and 
should be accounted for if they are found to exist. 

Summary of Ignition Transient Analysis 

During the course of this investigation an internal ballistics model 
of the ignition transient was developed using the Z-N burning rate model 
and incorporating the thermoelastic analysis. This was done so as to 
better represent the ignition transient than is possible by assigning a 
pure exponential pressure loading. Although the thermoelastic effect, 
which has now been shown to be small, is not augmented significantly by 
such treatment, the Z-N analysis may prove of some utility. Generally, 
the method is an adaptation of the ignition analysis technique presented 
in Refs. 6 and 27. but with the Z-N burnirg rate model ( .. 1, Initially, 
the thermoelastic coupling) included. The major work on the analysis was 
performed as a separate effort as a graduate thesis under guidance of 
this project and details will be presented in Ref. 8. Only a brief de- 
scription and illustration of the results is given here. 

The determination was made early in the analysis that the thermoelastic 
coupling has much less effect on the ignition transient prediction than 
does the unsteady Z-N burning rate. And since the thermoelastic coupling 
analysis needlessly adds to the computer time and storage requirements, it 
is not considered in the analysis of Ref. 8. The analysis compares pre- 
dicted ignition transients incorporating the unsteady Z-N burning rate 
law to predictions using the classical Vleille burning rate law (r « aP”) 
and to test data for the Titan lllC and Space Shuttle SRMs. The propel- 
lant properties necessary for the Vleille ignition transient are well 
known, but the propellant constants in the Arrhenius burning rate law, 
necessary for the Z-N formulation, are not so well known. Also, the 
flame-spreading speed is uncertain. Following the procedure of Ref. 8, 
these unknown quantities are adjusted to establish a good fit of the Z-N 
ignition thrust transient to the Titan IllC test data. The flame-sprcat.l 'p 
speed is adjusted so that the maximum rate of thrust rise for the V.i> x 
ignition transient prediction matches the maximum rate exhibi’ ’ '’V ; 

Titan me test data. This value of flame-spreading speed (4.iDO /se^i 

is then used throughout all analyses. The Arrhenius cons.;ri;*-t5 a :.>•* 

Justed to produce a good fit of peak thrust of the Titan ;\’C 



6 « 


transient test data. In the absence of any better Information, the same 
values of Arrhenius pre-exponential factor (740 in/sec), act vation energy 
(4900 cal/gmole), and flame-spreading speed are a^ed in the ignition 
transient prediction of the Space Shuttle SRM. 

Using these propellant constants, the Z-N and Vieille Ignition 
transient predictions are calculated for the Space Shuttle SRM and com- 
pared in Fig. V-11 to test data from Ref. 19. It is immediately perceived 
that both predictions exhibit a much smaller rate of pressure rise than 
the test data, thus underestimating the pressure throughout the transient. 
However, the Z-N model does exhibit a slightly higher rate of pressure 
rise than the Vieille model, and more closely predicts the pressure during 
the latter portion of the ignition transien .. 

It might be suspected from examination of the rate of pressure rise 
in Fig. V-11 that the flame-spreading speed is not of sufficient magnitude. 

In fact, MSFC, using a program developed by Caveny (Ref. 28), has analytic- 
ally predicted a flame-spreading speed in excess of 2’’,000 In/sec for the 
Space Shuttle SRM (Ref. 29). But this analysis also gave a predicted rate 
of thrust rise which is approximately twice that exhibited by the test data. 
This might indicate that the actual flame-spreading speed is approximately 
half that predicted. The results of using a flame-spreading speed of 
12,000 in/sec in the present Z-N analysis are given in Fig. V-12. The 
ignition transient pressure predictions are much Improved, and especially 
the Z-N model during the latter portion of the transient. Note that the 
pressure differential between the Z-N and Vieille models is greater 
with the increased flame-spreading speed. This is a result of the increased 
Z-N burning rates resulting from the more rapid rate of pressure rise with 
increased flame-spreading speed. 

Additionally, in Fig. V-12 it is noted that during the Initial 200 
milliseconds, the pressure is overestimated while the rate of pressure 
rise is underestimated. The overestimation of pressure might be at ’.riljuted 
to the fact that no delay associated with ignition Initiation is Incorporated 
in the analysis. And the underestimation of pressure rise may he caused 
by the assumption that the propellant grain is ignited adjacent to the 
igniter and progresses with time from that point. Bit, in fact, the pro- 
pellant Is probably ignited, after a delay, for some distance down the 
motor port, resulting In a higher-than-predlcted pressure, rise during this 
portion of the transient. It is planned to incorporate these cons Ida * Uious 
into the ana.lysls at a later time. 

Although there exists some uncertainty of flame-spreading speed and 
Arrhenius burning law constants, the Z-N theory does .show potential, as 
demonstrated in Refs. 7 and 8, for more accairately predicting the unsteady 
burning rates associated with Ignition and extinguishment transient.^. The 
Z-N theox'y has been demonstrated to be functional over a wide runge ol 
Arrhenius constants, but any additional quantitative conc.Uisions on the 
validity of the model must await more reliable values of the Arrhenius 
constants. 
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i-gure V-12, Effect, of increased flawc- spreading speed on the pre- 
dicted ignition transloiu for tl\e Space S' ittl.; SRM. 



VI. CONCLUDING REMARKS 


The continued research in Internal ballistic performance variation has 
produced considerable reduction In computer execution times for both the 
design analysis and Monte Carlo programs. Handling of Input data by the 
user has been simplified and the volume of printed material reduced. 

An alternative technique has been developed for predicting thrust 
Imbalance without recourse to the Monte Carlo program which is best 
applied to the talloff portion of the operation. A scaling relationship 
has been established to relate theoretical results to test results which 
may be applied to the alternative technique of predicting thrust imbalance 
or to the Monte Carlo evaluation. 

/d-^- tier -.1 evidence of the validity of the analysis of cirrular- 
perfor ;d ypellant deformation effects upon internal bau.listlc perform- 
ance h. bewii presented. The grain deformation analysis coupled with a 
new model of erosive burning seems to offer great potential for improving 
the predictability of SRM performance. The technique, however, needs to 
be further Investigated before the approach is incorporated Into the Monte 
Carlo program. Particular attention needs to be given to two-dimens .>nal 
nozzle flow and erosive burning threshold velocity effects. Although we 
have presented an approach for extending the grain deformation analysis to 
star-type grains, the method Is somewhat cumb . csome. The framework has 
been set, however, for a simplified approach. The more rigorous approach 
has also shown star grain Reformation would Se of small consequence fer 
many configurations, e.g. , the Space Shvttle. 

Finally, the extended investigation into the effect of strain rate on 
propallant burning rate leads to the conclusion that the thermoelastic 
effect -\s generally negligible for both steadily increasing pressure loads 
and oscillatory loads. No reason appears to exist for continued research 
in thermoelastic behavior. However, the capability to calculate the in- 
fluence of oscillating pressure loads on transient burning rates making 
use of the Zeldovich and NovozhJlov transient burning rate model may prove 
useful in combustion instability Iri/estlgatlons. 
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